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8UMMp .Y 


Tho irooults of Phaoo II flow ©xparlmantfj on tho opnn cyclo Has 
core concept are dlocuflaod In thin report. The flrnt pha8e<^> achieved 
flow configurations which gave high volume fractions for the inner gas, which 

“1' ono air as tho two gooS" both 
sj^ulatlng a propellant flowing around the outside and the fuel suspended in 
the center, volume fractions up to 33* in spherical configurations were 
obtained. Tho second phase of this work concentrated on improving the condl- 
tions whm gasos of different density wore employed, the central gas being 
no heavier. Again, both two-dimensional and three-dimensional tests were 
!! the two-dimensional oases, the round-shaped cavity was bounded by 
two flat walla on the vlwing ends. 

.. - 1 ^ largest tests were conducted on 3 ft diameter cavities, 

scaling up will be required to determine what the effects would be in the 

reactor which will have a cavity diameter of between 
of scaling from small models to larger sises was 
18 in. and 38 in. diameter cavities. The results indi- 
cate that neither the Reynolds nor the modified Proude number are satis- 
ractory scaling Indicles. The most nearly equivalent flow patterus for 

bMulexI wb«n . comproBi.. wm achUvad 

ijrnSJfiKJ" b ■Mtchlng Fraud, nuabar In one ca, «>d a match- 

xng Reyttold's number in another casee 

Several difficult and unusual aspects of conditions in the flowlne 
run iavestigated. For Instance, a dust injection test series was 

run. Sttt.h tests would simulate startup conditions in a cold reactor. Also, 
^slde down arrangements, with the exhaust nozzle facing up were tested, 
is puts the inertial forces in the opposite direction from those in an 

rocket. However, the upside down cases could be used in 
Jff essence, containment of the heavy central 

w ^®® ®*®®JJ®“*= these Upside down tests, with the principle problem being 
how to keep the central gaa ioray from the walls of the cavity. 

the method°« 4 f^J? normal dojmfiring directions, one of the main considerations is 
the method of dispersing the heavy central gas so that it occupies a large 
volume of high density within the cavity before it '‘falls out'' the exhaust 

11 ^^^® Sravltatlonal effect relative to the light 

propellant gas flowing around it. One mechanism that was found to be parti- 
e,. effective was a strong upward recirculation pattern for the outer gas. 
neS generated by the convergence of the annular gas streams 

outlet nozzle. Such recirculation was most noticeable In the three- 
duensional spherical configurations, and was undoubt><%dly dependent on the 
Shape of the lower part of the cavity. The upward recirculation tended to 
prevent the heavy center gas from dropping straight to the exit nozzle, and 
also dispersed the central gas. Since the central gas simulated uranium fuel, 
this dispersal r^^uced the self-shielding and enhanced the nuclear reactivity 
Of the system. This dependence of the shape of the lower portion of the cavity 
on containment that can be obtained with gases of different density is at vari- 
ance with results reported In Phase I of the study^^?. In those eases, princi- 
pally two-dimensional studies were performed and mostly on gases of the same 
density, where the gravitational effect was not of concern. 
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Pi«ur« 2 I*"® SE?” eo»axlal flowing gna cor© concepfe^^) g^own in 

*’«en studied oxparlroentnlly In a number of 
Inbor^orlea. Pattorna with high ratios of outer to Inner gas flo»^ ratoa 
studied in cyllndrlcal-goomotry configurations (3, 4). sub- 
whiEE i concentrated on spherically shaped configurations, 

closely resemble conditions which would be achlsvsd in actual 

'^^ich gave the moat clearly reactive 
“?«!«?« . .xporlmont. to study the roacto. phjoltrcta?. 

acteriatics Of a dilute gas core surrounded with hydrogen and a low abaorn- 
tlon moderator (heavy water) have also been reportadCsf??. since the ideal 

spherical, and since a gaseous fuel is to be used, 
considera;lons bear vary heavily on the types of flow patterns which 
Will be acceptable. In particular, the flow patterns, must be capable of 

operating pressure, which is 

S ? ^^®«® flow tests, con- 

ducted at Aerojet Nuclear, have concentrated on obtaining flow oonflguratlons 

pJssJbU^’SJth^e saa. Simulating the fuel, into as JaH 

tlons undM The goals for the condi.- 

wons under which the system will be operated are: minimum cavity pressures 

maximum propellant-. to nuclear fuel floS rate ratios, mlnSL ?eacLrSee 

apoUcILInf^M^i: *^®^®® *^® corresjond to the low thJuS’ 

fEturr^ ®*^® Tor nuclear rr :kets of the not too distant 


\ 
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Fig* 2*1 Open Cycle Gas Core Concept 


3.0 


TS3T APPARATUS 


The fXow teat npporntua waa dealgncd ao that modificationa could 
he made enaily to the cavity ahape, to Ita walla, to the lolTJ^ lZuZtt 
finM* Injoctora, and to the oKhauat noRBlo. The Important conald- 

dlanaaombly and roaoBambly, Fine tuulnR of the flow 
dlatrlbutlona into tho cavity and muaauromonts of the flow parametora 
were of secondary importanco. However, tho total flow rotos for both the 
inner and outer gaaea were moasurod with hiHl» accuracy. The boalc apoaratus 

blower wltlia rated ca;acltj ori^OO^^cfn^- 
vidod the main outside gas flow. Air was used as tho outer gas In all caSes 

??«ur thi IsTnVi ® ^“"80 teat bL. In the ’ 

if Jr iu f diameter spherical cavity can be scon. This waa constructed. 

in t. « Phase l^reoorfc^^^ and is the identical sphere that was described 

bL « report'- K %e central gas flow was provided from a smoke mixing 

a 2’ireLr?aSo thence Ih^J^gh 

° ^ central gas injector inside the cavity. Tho mixing box 

which wiuis^ volume, and hence provided a relatively large accumulator 

con..^ant smoke density for reasonable periods of time. The 
ohitoielic?rfi® ^ox near its top, was continuously monitored by a 

dewlt? f« all^'Jiiti was used to establish a uniform smoke 

cha^bw ^ orific! the smoke entering the pipe on its way to the teat 

amber. Orifice flow measuring sections were Installed in both the 8 in and 

i!«f!n ? ^ versatile, allowing cavity to be 

a dowflrlng or uttflring eoiU*lgurations,*and provided 

®cre even distribution of flow into the test cavity region. The large 36 
in. diameter plastic sphere is shown in Figure 3.3. The dimensions of^thls 

sS«e “‘i ? ” “P ftom th... of the U 

aS c«rh/rfr a «iaot differences between the two spheres. 

Ma« sphere did not have a circular croaa sectlou but 

facilitate conatructlon of such a large 

fM perfornl^g”?h°^m^?5?® ^ diameter rings waa difficult, and facilities 
tor performing the molding were generally unavailable. Furthermore, several 

off foi'’^hr?e8?rwhfh‘* fsf! ®P^®’"®* ^^^®®® blocked 

Ing drLliw of^hf if if f ® up conditions. For comparison, engineer- ■ 
3?4b. * of the 18 In. and 36 In. spheres are shown in Figures 3.4a and 


4 ..u 44 ^®*^ *^^® ®u**ualled "two-dimensional" tests, the cavity ahaoe was flat 

facilitated viewing attd interpretation of the tests and made the construction 
of a number of different cavity shapes relatively simple. The flat end walls 

somewhat. However, this end perturbation extended 
y a few inches from the viewing wall surface, and at least qualitatively did 


SR)*^cf Itlef tested incorporated downward firing (downward exhaust 
riott f t f 4 gravity to simulate the effect of forward accolera- 
fofa engine in space. Upward firing can be utlllaad in practice 

f i,\ **®''®*^ ^ application, or possibly for simulating 
rocket with a thrust reverser unit. 
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not c6wppoini3e the interpretation of the test resuite, Clooe-up photographo 
of the two contignratione uoed in the two'^dimenaionai tenting are ahown in 

36 in. cavity). Again, each 
«‘«ct acali iBodol of the other. Each containa an outer 
cavity wall with perforations, and on this wall were attached 24 perforated 

directing the flow in an tangential direction along the cavity wall, 
iguro 3t7 ahowQ croai^sectloftal details of the 36 in# test cavity for a -typical 
sanL^J®!; T discovered in Phase X of this tooting (Reference 1, 

J J . lu tranoplratlon cooling holes in the louvers should be 

kept to the minimum necessary to provide adequate cooling. When this la 
done, the transpiration flow is small and has a nogllgiblo effect on the 
Thsreforo, for many of the testa the tranoplra?5on SooUng 
''®f® the tliree-dimonsionel transparent 

the flSi.r transpiration cooling holes, only louvers to direct 


!®“ hy a blower, and sinoe pressure 

®PP®tatU8 ware fflinlmal, choke flow conditions could not 
i® tolerated. The nossle areas were at least 1 square- inch. Phase I teat- 

had little or no effect on the flow patterns established inside the cavity. 
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Fig. 3.4 (b) Vertical cross section oi large, 36-ln. diameter •’spherical” cavity 

10 











12 





Cross section of 36 In. two •dimensional test cavity used for 
typical test series. 


4.0 


TliST PROCEDURES 


ThP procedurea fot operating these tests capltailaed on the flexibility 
and ease of expcilmentol changes. The front face of sJther test box could bo 
easily removed to make alterations on the cavity arrangement or to the inlet 
injector. After the front was replaced on the box and sealed, constant flow 
was then oatabllshed for both the outer and inner gas circuits. The flow rates 
wore monitored visually on differential pressure motors. A standard series of 
even numbered flow rates for the outside gas were used for most of the test. 

These ranged from 50 cfm to 600 efra on the 18 in. cavities to 100 to 1100 cfm 
on the 36 in. cavities. The inner gas flow rates were established so as to give 
an even numbered ratio with the outside gas. For the initial experiments, these 
flow rate ratios were even in terms of volume flow rates. But for the latter 
phases of the experiments, even numbers were chosen for the mass flow rate ratios. 
Flows were regulated by means of valves at each gas source, i.e. , the blower for 
the outer gas and the pressure regulator for the inner gas. The only adjustments 
provided beyond this point were at the walls of the cavity, where Scott foam or 
plastic sheet was used to partially or fully restrict air flow through certain 
portions of the cavity wall. The basic data for each test were the two flow 
rates and the photographs that were taken of tlie visible flow patterns. 


Among the types of supplementary data obtained were flow velocity 
measurements within the cavity. These have been conducted to date only on the 
two-dimensional configurations, because of the difficulty of obtaining access 
to the interior of the three-dimensional ca'^ities. These velocity measurements 
were made with hot wire anemometers. In general, the gas velocities were inde- 
pendent of the inner gas flow rate. However, in the case of very heavy gases 

gravitational effect perturbed the main flow patterns 
centerline of the cavity. In such cases, flow measurwnents were taken 
with the hot wire anemometer with but a single gas flowing, since gas mixtures 
will perturb the output of the anemometer. However, to determine if the heavy 
did todeed have an effect on the flow patterns, smoke traces were 
inserted along the centerline and in other sensitive regions of the cavity to 
deteralne the exact direction of the flow in these areas. (The anemometer 
itself would not determine the direction of the flow} it was only capable of 
determining the line of the maximum velocity or gradient.) Other, and perhaps 
more Important, supplemental information was obtained by scanning the photo- 
graphs with a densitometer in order to determine smoke densities, and hence 
densities of the central gas, throughout the cavity. To obtain satisfactory 
pictures for the densitometer scanning, time exposures were taken. These 
averaged-out the time dependent fluctuations in the smoke density. These time 
wposure photographs were exposed such that the linear portion of the log 
density exposure curve of the film was utilised. Selection of the proper 
exposure times ore discussed in Reference 1, Section 4.0. 


%e time exposure photographs were scanned on a Joyce automatic scan- 

obtained radially across the cavity at a number 
of different longitudinal locations. The resultant traces were normalized for 
smoke densities varying between 0 and 100%. The 100% value was established as 
the density at the outlet of the Injection nozzle. For the two-dimensional 
cavity tests, the effective density through the cavity could easily be inter- 
preted directly from the smoke density traces, and the isodenelty lines could 
be directly plotted accordingly. 
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In the case of the three-diraenalonal cavity, with chord lenatha 

interpretations could not be directly 
Thf unfolded from the observed density traced 

equivalent to solving for the coefficients in a 
foJ rL "»ethod used and computer program 

k seherao are shown in Appendix C of Referonce-1. Further- 

reduction program so as 

the varying attenuation of the walls of the PlLiglas sphere 
rafuiLd^A^r^ik”® J" ijshting conditions. Thus, tho three-dimensional tests 
Po? involved process to obtain the radial smoke density function. 

leLJh flid ^^® ''s® simple, since chord-viewing- 

if background effects were essentially constant. However, the e- i wall 

J®** ^*»«/wo-dimenslonal results and made these density measure- 

simplifications were made in manrcLes 
in the two-dimensional density scans so as to obtain a larger quantity of data 

aL configurations, a single vertical 

volu^A ^5®®? ''®® estimates of the integrated 

8carcaferii?h comol»J^®f selected samples of such sLplified 

percen? a? ?hfmos?! scans showed that the errors were only a few 

accuracy of the densitometer scans is quite good, far 
of the °^ J®P*^°<^J°ihillty. The estimated accuracy on the integration 

tliiJv ?rJi estimated to be about ±10%. The principle uncer- 

cJnditlM of inS®^i°f®®!' however, is establishing a base line 

condition of 100% smoke density. Since this could only be determined by a 

measurement directly at the surface of the injection noazle, the assumption 
smoke density at this location was indeed 100%. Gas samJlL 
ImL location with the configuration running on air and 

thit ®“°^!* results of these gas samples indicate 

Th^^arr^A? A«*^/® ®^^^ ®? surface of the injection nozzle is nearly 100%. 
dlffere^i indicated 87% density, but a portion of this 

the JolS^e I attributed to the sampling technique, which sucked 

o ? ^ l! J sampled in rapidly. The net effect of all uncer- 

to^SltSn^20rof IhT estimated the volume fraction results are accurate 
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!» . 0 ;n4ST JLliHUJVj'H A ND UIKCUtiSl QN 

mn.jor Hoctlon of tho report doais nnparfttoly with tho varloun 

proRram. Each of tha major napacta of the variableo 
^ttLnt^hn! conflRuratlon arc dealt with in tho aubncutlona. L ^ 

attempt haa boon made to nrranR® theao rcnultfl Into a npoclfle order; tltev 
are liatod In the followiiiR order Cor convenience on3y, ^ 

Inlectlon Nogzlo Effoeta 

fflr ^ thcHQ Btudlea^^^ clearly showed that the Injection noazle 

V alternatively, the method of injectinij the center gao. 

JL ^ acceptable flow patterns and densities for the InSor 

gas. It should not be assumed that there will bo but a single Inieotlon method 

ake for simplicity of design of the actual reactor, further studies will be 
necessary before the final configuration can bo selected. 

«<.. ^ assumed that In such a roc'.ct, either pel.l,eta of metal or a thin 

uJnn «jaalum metal will be Injected rapidly into the cavity -*rom the top. 

lUerauf^DlJde 171 Interior of the cavity, the metal will 

this newlv formed ®°“® mechanism must be present to disperse 

in creaSL n ®? regions of the cavity where it Is most effective 

nuclear chained reaction. Thus, the non-nuclear flow testlim 

very^lL introduce the center gas into the cavity with 

very low velocities. This essentially takes little or no credit for any exolo- 

velocity which will result when the metal comes in contact wlth^the 

deslSd LnnLoL ^ flowing outer gas, is one method of achieving the 

° ®i®P®rsal. However, too much entrainment of the center gas Into the 

Jjfwairj ^r Sas being quieSy ^ep? abound 

involves settina SfnL The second cSnaideratlon 

j ^ ® flow patterns of a non**turbulent nature In the central 

regions of the cavity, flow patterns which will tend to buoy up the heavier 

under^th^^< buoyancy forces, a heavy center gaa tends to fall 

gravity directly to the exhauat nozzle, resisting what- 
ever other dispersal forces might be present. 

5*1*1 I'Two-Dlm ensional'* Ini ec tor PesienB 

inlector ^^‘‘Tigurations represent an unusual problem in 

meieli be a basic configuration was intended to 

loif Jlewinu iord w spherical arrangement. However, the 

thrLahout the two viewing plates should represent 

InWlornoaaiP^fi^ ? uniform flow pattern. This is not possible unless the 
tSo JlwlnrsuiLoi! ^ uuml^lly equivalent to that between the 

•nTn<l.« «L 4 ^f ® Tong injector, plus the desire for 

minima, and virtually zero, injection velocity are generally incomoatlbla 

fiSit! ?h! ^ ? f®?* ^ pressure drop on the other hand will result in 

faa^^Ld«^r velocities. Conversely, if no pressure drop is employed, the 

irthrend of“tho ®J>'®Tght dojm the pipe until it reaches an obstruction (sJch 
end of the pipe) , at which time the gas will come through the Inlectlon 
no«l. op.„l„g. into the cavity. Figure 5.1 .hew. the geSS'JyjS cF 
be expected under those circumstances. If flow is introduced frL both ends 
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patterns aiwUar to that shown in part B of Figure 


<no K^oci^**”*** improvement in the above patterns can be obtained by fashion-" 

the flow an desired inside the injector 
noaslos. Some typical plug designs are shewn in Figure S.2. These generally 
achieved the type of uniformity in injection pattern doslrod. Their one dis- 
advantage was t^t when designed to work on one typo-ef injection gas, the 
nozzle did not function os desired on a different density gas. 

i-K *.4 “bove described difficulties in achieving uniform injection over 

length of the injection nozzle, plus the oompronises already inherent 
shoJld configurations, dictates that the two-dUnenaion^l resultJ 

?o In a qualitative manner. The two-dimensional configu- 

ration is much easier to interpret and to understand in flow patterns, visual* 

fied chan, is the case with the spherical models. A further advantage is relative 

that tL'^ti^diSeMLnai^fillf/'^® louvers can be varied. It is for these reasons 
cnat tne two-d^ensional configurations were employed, and developed into a basic 
arrangement which gave the most satisfactory flow resiles. 

Spheric al Inlector Dealena 

Iona and short cylinders, approximately W/2 to 2 inches 

coSligSjatlinf f injection nozzlea in the spherical 

contxguratlons. However, these also had similar disadvantages to the oloe con- 

o^Z^^lZ two-dimensional cavities. The flSJ^Jended to lilt 

the pipe. Again, the Insertion of especially made 
or baffles was able to obtain a more uniform distribution of flow over 

ini sc'h on *^^® Injection nozzle. However, a more aultable pattern for 

S u*f^? arrangement would be a spherical Injector containing perfora- 
are shown In Figure 5.3. The 3 in. diameter injector 
SJ ?n^ the 6 in. diameter injector waa used in 

awe^biS* r spherical injectors dls- 

aas^bled. The Inner perforated cylinder formed the pressure wall which 

helped create relatively uniform dlsperaal over the spherical surface.. The 
second or outside wall then allowed the gas to issue forth into the cavity 

IqTII injection velocity. The outside walls of the spherical injectors 
®°nei®;ins oi nominally 0.125 in. diameter holes. These 
ttjec tore generally worked quite satisfactorily when air waa the injection 

fM^the epl'e^^icnl surface.^ However, 

*®f®®; tendency still existed for the gas to 
Selief ^ portion of the sphere. Since a uranium 

^®i»8 injected ir.t the cavity would also have a 
prefer^tially velocity towards the nozzle from both its initial velocity and 

be°r6alistlc?^^*^**^^”* effect, these effects with the injector are believed to 
5*1*3 Effect of Inlector Position 

^nlo^ « position In the cavity will determine how mubh of the 

inlet gas can be pulled out by velocity shear forces originating from the 
fast moving outer gas flowing along the walls. The original results of such 
investigations are given in Section 5.1 of Reference 1.* In brief, these 
tests showed that an optlmifio position for the injection nozzle can be found 
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From inlet air 
and smoke source 


From Inlet alr-^nd smoke source 




Top picture shows effect when smoke was let into injector through single oloo and 
bottom picture 8how8 effect of double Inlets Into the injector! I JjStor cSnllst^^ 
of a perforated metal cylinder covered with a layer of Scott foam. %e end of^the^slngle 
inlet Injector was a solid disk. (Note, rectangular box outline has no significance.) 

Pig. 5*1 Injectors used for Two-Dimensional Smoke Teats 
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Spherical injector details 



which glvcR maximum volume fraction of the central gas in the cavity. When 

Involved for the central gaa, It becoAea 
m.. K provide aufflclent dlaperaal forcoe without alao entraining ao 

rignlflcantly**”^*^”^ outaldo gaa stream that volume fraction suffers 

tests In the spherical arrangement involving hoavlor than air 
gases for the center gas, an unexpected buoyancy force resulted from the re- 
circulation patterns produced in the Interior of the sphere. These recircula- 
tion patterns are shown in Figure 5.5. The effect of these recirculation 
patterns on the central gas is shown in Figures 5.6 and 5.7, for the 18-inuh 
and 36- inch spheres, respectively. The upward recirculation with air as the 
center gas is quite evident, and even with argon there is evidence of signifi- 
cant buoyancy, enough to even cause the argon to circulate upwards. As can be 
era, the position of the injection nozsle is not nearly as critical with this 
spherical arrangement as it was in the cylinder core arrangement. The latter 
did not have nearly as strong as recirculation pattern and hence needed to rely 
more heavily on the velocity shear effects. 

The net effect of such recirculation patterns is not only additional 
buoyancy but additional mixing of the inner and outer gases. However, mixing 
of hydrogen with the uranium gas is not necessarily a deleterious effect, unless 

uranium in the cavity. What is Important is to 
attain the highest possible mass of the central gas distributed in as large as 
possible radius inside the cavity, in the constraints of the desired low flow, 
rates of outer tubular gas. 

Efiecfc of Preferent ial Direction of Center Gas imectlon 

, 41 . Though it has generally been assumed that injection of the center gas 
will occur from the top of the cavity, this need not necessarily be the case. 

In fact, since the pellets or fine wire will be injected with a rather high 
velocity, directing this velocity toward the exhaust nozzle may be considered 
an inappropriate design. Therefore, some brief studies were done on injection 
center gas from other positions along the cavity wall. In particular, 
upward injection from a mechanism located near the exhaust end of the nozzle 
was ^amlned. Figure 5.8 shows comparisons between preferential upward rad 
preferential downward injection. 

As can be seen, upward injection offers an additional control mechanism 

acceleration forces which tend to move the heavy gas 
flexibility would not exist in the real operating gas 
core rocket. There, no mechanisms could be positioned inside the cavity, at 
4 4 ®^.. 4 ®®^ ^haust nozzle end. However, it might be conceivable that the 

injection velocity of the pellets or of the fine wire could be varied signifi- 
®o that the effective injector position (i.e., the point at which 
in id ^^h*^ different longitudinal positions 
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Pig, 5.5 Typical cavity recirculation patterns 
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in. 3-D, argon, injector at 2 in., 5 in., and 10 In. do»n 
(Some louvers blocked with Scott foam) 

18 in. spherical test cases show^g recirculation effects arising from 
varying 3 in. Injector position 


















Downward Injection of amoke (Sam® flow conditions) 
Fig. 5.8 Effect of preferential direction of Injection 





5.2 


Scollng Studlea 


^ problem of primary Importance In developing analysis technlaues 
tiLUT^ technology la that of acallng from small test conf^gSra- 

cavitlofl. In an effort to determine the applicable 
d£Seiai«i«r cavity, tests were conducted on th^two- 

7 ® oavltlM war. of OBaantlally Idantlcol geonatry oven 

2* f*” ‘nlnnlnn 'Itoonslon.. The 36 la. apharlcal cavity wa"a factor 

“f “ circular croas aoctlon as for the wall one. 

should “'■”‘”•"*“1 '““1“ 

f«». kaad-u 5*^ bo 5a 12 show examples of some of the scale-up comparisons 

iS «irS^”'‘rSl' " J?' "T' *"• <“* Saa was aJr la 2Il“«22?) 

eJnaf-Ht 1 ^own, the equation of motion in fluid flow Is dimensionless for 

well ^[oude and Reynold's numbers. But as Is equally 

la ^ / “®v® 4T to maintain dynamic similarity If the^slze 

%ere?Sra ’t!:; characteristics can be manipulated also! 

R^JJ^ld^rn^w ® comparison example In which the 

I” ^*C'™I“I”8 the scaling for these numbera. It was aasuned that the 
l^t velocities vary Inversely with surfwje are. chaige. Hence. In scaling 

Se2 L‘22i2t22t ?2r /I“«“ Jlmenelons doubling, the HeynoldU 

aT 5 r * ?i5"™iJ2r:L°:-in^for^J2S 

oen an a. ^ scett, neither the Reynold's nor modified Froude numbers 

SdJSte\«® aatlsfactory scaling IndicS. However, the res^Us geJI^^SL 

^wde an^RevS^dU similarity somewhere between the matching modified 

Sf! and Reynold's number conditions, but favoring somewhat the Reynold's 

Setter c!J!ZT^ “P ft diameter tHV 

to 8 or 10 ft Ai f losical manner and that similar extrapolations 

to B or 10 ft diameter cavities can be made, but with care. 


5.3 


Buoyancy Effects 


ai>olleAtiiJ«^!°Jj®”J hydraulic parameter In rocket engine 

dMiioJ It af'her. This dimensionless number Is empirically 

orjjf fSel!n?Le?w acceleration In ouJer ^pace) 

on cne tuei propellant In the engine chamber. In the aas core reari-nw w»r.iewl 

engine, the replar Froude number has to be modified to account for the differ- 

f * propellant densities. ^tri^enu^LteJ Sat the 

o If led Froude number, called the buoyancy number, be given by 


B 


/ 


Bk 

2 

V 


h£. 

Py 


27 



upper pictures, 18 In. cavity, 100 cftn/0.72 cfm and 200/1.44 
Lower pictures, 36 In., cavity, 720 cfm/2.88 cfm and 400/2.88 

Pig. 5.9 Scaling series with two-dimensional cavities 
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Upper pictures, 18 in., cavity, 100 cfm/1 cfm and 200/2 
Lower pictures, 36 in., cavity, 720 cfm/6 cftn and 400/4 

Fig. 5.10 Scaling series with two-dimensional cavities 
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upper pictures, 18 In,, cavity, 100 cfm/1 cfm and 200/2 
Lower pictures, 36 In. cavity, 600 cfm/6 and 400/4 

Fig. 5.11 Scaling series with spherical cavities 
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which Ifl the square root of the inverse sf the ususi Proude number, times 
a factor for the relative density effect, Expefiroentation has shown that 
as the relative densities increase, fuel containment decreaaes, roaklns the 
buoyancy number proportional to the mass ratio of the inner and outer gases, 
and inversely proportional to the fuel containment. 

The buoyancy numbers wore calculated for a large sorion of 2-D tent 
CMOS for both the 18 and 36 in. cavition. The b in the formula was taken 
in all CO0O8 to be tho nominal distance from the injoetor to tlw exit noealo. 
a reasonably roprosontatlvo measure of tho distance a fuel particle travels 
within the cavity. The velocity fiorm In the formula Js squared, which fact, 
coupled with poor correlation hoped for between buoyancy and fuel contain- 
choice of "combined" buoyonoy numbers, using tho prodtict 
o two differently determined velocities for tho velocity squared term. Those 
calculations ate sunnnatlzod In Tables 5el and 5*2e 

Initially it was assumed that tho velority to b« used In dotorminlng 
buoyancy should be the injection velocity of the inner gas. The measure of 
ruei containment to which the buoyancy numbers should correspond is taken to 
K« 4 fraction, the ratio of volume of the cavity that would 

'>^®’f® collected into a region of 100!!! 
cavity volume. Figure 5.13 ahows this initial buoyancy 
ttUffloer, called the Injector Buoyancy number, vs the fuel volume fraction. 

jmiie a general trend can be seen, better correlation was hoped for, so additional 
buoyancy numbern were postulated. aaaicxonai 

Figure 5.14 shows the Recirculation Buoyancy number, which uses 1/10 
of the outer gas flow velocity for the velocity squared term. *^18 scheme 

buoiiSS fraction of the outer gas reclrcul utes within the 

cavity, buoying up the Inner gas. Correlation is poor In this case also. 

,iai.. 4 ..^u depicts the results of the Combination Buoyancy numbers, 

by combining the injection velocity times the recirculation velocity 

III ^^®^ocity squared term. Again, no distinct correlation between buoyancy 
ana volume fraction was observed. ^ 

no definite correlations were found with the previously mentioned 
hJ*?kJ ® buoyancy number, an Area Buoyancy number was defined 

JL velocity to be used by assuming that the volume of the Inner 

time crosses the m iplane turbulence area of the cevity. 

Attain® a ^U®^* ^^®f® displayed in Figure 5.16. 

point gM, and 


The last buoyancy number shown, the Reaidence Buoyancy number. Figure 
5.17, n«ed a velocity derived from the distance from the Injector to the exit 
nottxle divided by m average fuel particle residence time. The residence time 
was determined by dividing the volume occupied by fuel by the volume flow tatee 
to give a reppeentetlve time that a fuel particle recirculates In the cavity. 
Here, as In the preceding oases, it can be seen that, in general, as buoyancy 

c^ld*^be**idMtifi fraction decreases, but no optimum buoyancy numbers 
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Fxg. 5.13 Cavity volume fraction vs injector buoyancy number 










Fig. 5.14 Cavity volume fraction vs recirculation buoyancy number 









Fxg* 5«16 Cavity volume fraction vs area buoyancy 'number 









Fig. 5.17 Cavity volume fraction vs residence buoyancy nunber 
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Dup tP thp length of time needed for reduction of aphericnl acen 
casee to volume , actione, little aphericel buoyancy data were amaaaod. 
Becauee - 0 £ the euojeotivity of the volume fraction determinationa from 
the deneitomoter aeon caaoa uaed, the data in thia aection ahould be» token 
to be qualitative only. 


S.4 


aflAocity Meanuromontn 


Velocity moaeuromontR slmilor to those reported in Phase I of this 
Study woiTQ tdken in tho 3 ft dloaotot 2-D cavity using a hot wire anemomstote 
Ueing outer propellant flows of 200, 400, and 720 cfa, tho cavity was mapped 
for both velocity gradient and direction,* However, oven at 720 cfm outer 
flow conditions, the sign of the grodient would not be easily established 

in!: cavity. In that interior region, velocitlea 

were less than 25 ft/mln. At these low velocities, tho usual direction indi- 
cator (fine thread) did not respond. Later attempts using smoke tracers were 
successful in showing that the recirculation pattern inside the cavity was 
osttTOSibly upward, away from the exhaust nozale, all the way out to 2/3 of the 

-.u Since hot wire anemometers are perturbed by different density gases. 

mostly for air as both outer and inner gas. Due to the 
Accessibility of the Interior of the three-dimensional spheres, no velocity 
measurements were made for these cavities. 


5.5 


Sweep-Out Rates 


Most of the testing repotted thus far has been concerned largely with 
p pcAetcfs of- the central gas region under steady-state flow- conditions 
tor both gases. If there were no sweep-out of inner gas by the outer gas, a 
very low fuel injection rate would be needed to replace fuel lost to the 
system through fissioning and neutron capture. Therefore, very low, or even 

tSts c^ducted* primary Interest in the sweep-out rate 

The recirculating center ges region exists, whether there is enough 
tracer smoke to render it visible or not. The size and character of this 
central region depends on the cavity and flow parameters for the particular 
coi^Aratiott under study. It seems plausible, however, that there should 
exist a minimum effect gas/lnjection rate below which the central turbulence 

change. One way to determine that minimum effect 
central gas once steady state has been reached and measure 
u Which the emoke Is swept item the cavity » this 8wesp«-out tate 
should be Aependent of the central gas injection rate at values below the 
mittlmutt effect rate, and should be a measure of the minimum fuel loss rate 
for a particular configuration. 


Figures 5.16 to 5.23 show the maps obtained from these measurements. The 
outer fast moving layer can be seen to be of relatively constant thickness, 
the velocity peaking at around 14 in. radius. 
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Fig. 5.22 Typical velocity profile map of 36 In. two-dimensional test cavity 
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For the sweep-out rate tests » air was used for both inner and outer 
gs8» with the injector 10 in. below the top of the 36 in. 2-D cavity. A 
photograph was taken at steady state operating conditions » and the inner 
gas was shut-off. A timed sequence of photographs were then taken until no 
>f®t>alned in the cavity. The negatives were scanned across 
midplane with the densitometer, and an "eye-ball" estimate of the average 
smoke density for each picture in the series made. These average densities wore 
plotted vs decay time as shown in Figure 5.24. Average sweep-out rates were 
determined from these plots and graphed against initial center gas injection 
rate in Figure 5.25 and against outer gas flow rate in Figure 5.26. 

u It* sketchy and qualitative nature of these results do not appear to 
e sufficient grounds for making any firm conclusions* The results were not 
totally self-consistent,® the sweep-out rates were not consistent, let alone 
constant, during decay time. No minimum effect rate was determined, and If 
one exists. It Is at a higher value than any tested. The last two figures 
indicate that the sweep-out rate seems to be a rather weak function of initial 
center gas Injection rate, and a stronger function of outer gas flow rate. 

No noticeable change In the size of the central turbulence region was detected 
In the course of these tests. More exact Information on central gas densities 
^st be obtained before additional conclusions can be drawn from this type of 

t»d6 1* e 


Upflrlng Tests 

The results of running the test cavities upside down, so that the 
effect of gravity Is reversed arc shown in Figures 5.27 to 5.29. Good 
containment of all Inner gases was dmaonstrated. At low outer flow rates. 

It becomes difficult to prevent the heavier Inner gases from settling down 
very near to the cavity walls, which would create a problem in an operating 
^ctor because of the heat. A configuration of this type could be used In 
MHD or heat source related applications, while the problems Involved In pro- 
pellant reversal make this configuration Impractical for rocket engine appli- 
cations. 


Dust Injection Testa 


The effects involved in solid particle, or dust. Injection methods 
were Investigated. Such schemes would be of utility In start-up of an operable 

effect of upward vs downward initial scattering of the 
dust la ahowti In Figure 5* 30* A small variable speed rotating propeller fan 
was the scattering method used for the dust. A small amount of experimentation 
was done dust composition, that Is flour, bicarbonate of soda, AI 2 O,, and 
talc were all tried# Talc seemed to give the best results# 


One problima Identified with these dust injection schemes can be 
identified from the background pictures shown In Figure 5.30. That is the 
dust tends to settle out onto the walls of the cavity. More effort will be 


® Pointing out once again the very complex nature of the Interactions within 
the flowing gas cavity. 
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4 . 


Center Gan Region Smoke Decay 



Peripheral Gas Flow ■ 130 cfm 
Initial Center Gas Flow 

Slope (SOR) 

• 4.0 c£a 0.0570 

A 2.0 c£m 0.0585 

^ 1.0 cfm 0.0550 



Decay time (aec> 

5.24 Cavity smoke density vs decay time 
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Sweep-Out Rate (Ap/At) 
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Pig. 5.25 Cavity sweepout rate vs Initial center gas flow rate 
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t Bate (Ap/At) 





Fig* 5.27 Comparison of upt^ard exhausting cavity case to dovnvard 
exhausting cavity 





5.28 Camparison of upward exhausting cavity case to downward A^haiigr-tTig cavity 














Blank 


Upward scattering rotor 



Fig. 5.30 Oust Injection series showing effect of scattering dust both 
upward and downward 
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Id iWL'rt tirntna 


ATR-TALCUM POWDER 


«v3^i) ' 




Run i?"9"71, 3 300, 600r.|0 

Injection no?.?jp one 
Inch from top, Two- 
V Inch exit no?,«jp with 

S\ no other fittenuntiori. 

' Contour interval In 
' \\ % of center powder 

\\\ i!oncentration, reliitlv 

^ \ \ 100^ at l.jjeciion 

\ \ \ noa?.Io, 








Iso-density contours in spherical cavity. 
Integrated center powder concentration 
■ 0,l6 of cavity volume. Outer flow 
600 cfm, inner flow auger setting on 10. 


Isodenslty contour msp for typical dust (talcum powder) 
Injection test 



Fig. 5.32 CatR()uter Solution Model I 

3SX Volume Fraction, 0.93 Radius Ratio 
Reff ■ 82.35 cm. , Total Fuel Mass • 119 kfl. 
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Fig. 5.33 Computer Solution Model IX 

35S Volume Fracitlon* 0.93 Radius Ratio 
Reff ■ 108.78 un. » Total Fuel Mass ■ 119 kg 





Table 5.3 


Model 

I 


IV 


VI 


IX 

X 


35% Volume Fraction Model Specifications 


Region 

No. 

Material 

No. 

No. 

of 

Mesh 

Points 

1 

1 

20 

2 

2 

5 

3 

6 

5 

4 

7 

10 

5 

8 

15 

6 

10 

25 

7 

100 

10 

8 

200 

10 

9A 

201 

20-40 

9B 

201 

40 

1 

7 

11 

2 

6 

2 

3 

2 

2 

4 

1 

7 

5 

2 

3 

6 

6 

3 

7 

7 

3 

8 

8 

5 

9 

10 

9 

10 

8 

10 

11 

10 

25 

12 

100 

10 

13 

200 

10 

14 

201 

40 

1 

10 

30 

2 

8 

15 

3 

7 

5 

4 

6 

5 

5 

2 

5 

6 

1 

20 

7 

100 

10 

8 

200 

10 

9 

201 

40 

1 

35 

80 

2 

100 

10 

3 

200 

10 

4 

201 

40 

1 

1 

50 

2 

100 

40 

3 

200 

10 

4 

201 

40 


Cavity % 
Volume 


Outer 

Radius 

Region 

Thickness 

Occupied by 
Rep, ion 

84.75 

84.75 

17.2 

92.28 

7.53 

5.0 

96.00 

3.72 

2.8 

104.02 

8.02 

6.8 

121.66 

17.66 

19.1 

145.04 

23.36 

35.3 

152. 4 

7.36 

13.8 

162.56 

10.16 

— 

198.12 

35.56 

— 

243.84 

81.28 

— 

56.00 

56.00 

5.0 

60.00 

4.00 

1.0 

65.00 

5.00 

1.8 

96.00 

31.00 

17.2 

100.00 

4.00 

3-. 2 

102.00 

2.00 

1.8 

104.00 

2,00 

1.8 

110.00 

6.00 

5.8 

120.00 

10.00 

11.2 

130.00 

10.00 

13.3 

145.04 

15.04 

24.1 

152.40 

7.36 

13.8 

162.56 

10.16 


243.84 

81.28 

— 

107.71 

107.71 

35.3 

124.41 

16.70 

19.1 

129.39 

4.98 

6.8 

13.33 

1.94 

2.8 

134.67 

3.34 

5.0 

145.04 

10.37 

17.2 

152.4 

7.36 

13.8 

162. ;;6 

10.16 

— 

243.84 

81.28 

— 

145.04 

145.04 

86.2 

152.40 

7.36 

13.8 

162.56 

10.16 


243.83 

81.28 

— 

107.40 

107.40 

35.2 

152.40 

45.00 

64.8 

162.56 

10.16 

— 

243.84 

81.28 
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Tablo St 4 


Material Spociflcatlona 


Material No. 

Comoosition 

AlmJDena lti.cn (x lO"^ atomn.'cm) 

1 

100% U 


1.958 

2 

90% U 


1.762 


10% H 


0.5905 

6 

50% U 


0.9789 


50% H 


2.953 

7 

40% U 


0.7832 


60% H 


3.568 

8 

30% U 


0.5874 


70% H 


4.134 

10 

10% U 


0.1958 


90% E 


5.315 

35 

40. U U 


0.7950 


59.4% H 


3.489 

100 

100% H 


21.090 

200 

100% BeO 


7280. 

201 

100% O 2 O 

+ 0.0022%H 

3310. 



Fig. 5.35 MONA. Computer Solutlona — K vs Radios Ratio 








Volume Fraction at Constant Radius Ratio 



Fig. 5.37 Computer Solutions - K vs Volume Fraction at Constant Radius Ratio 




necessary to perfect these schemes to protect the cavity wails during operation 
nf solid dust is to be injected during startup. Figure 5.31 is 

bsokground reading was 

fni S. frames exposed before and after the tests to account 

th« walls. Because the talc Is such a fine powder, 

air about 1.6 gms of talc per about 10 tpns of 

schcroo is*^to°operatc**^^^ pteasures will have to be investigated if this startup 


5.8 


Criticality Calcul ationa and Selectod Reaulte 


comnuter nuclear criticality to flowing gas parameters, 

diffusion and transport theory computer calculations were 
hrLnf!!i?^? ^ idealised base configuration was arrived at 

M J ! isodensity contour plots drawn in the flow testing. Varia- 

tions on this base were performed to determine qualitatively the magnitude of 
nuclear criticality changes associated with possible fluctuations in the central 

models are shown in Figures 5.32 to 5.34. Some 
listed in Table 5.3. The material 
frrm. Table 5.4. Ilydrogtin compositions used were obtained 

trom Fatcht ^ for 500 atm and the appropriate temperatures. Virtually all the 

®*^ *^*^® temperatures and pressures of interest in the 

densities were obtained from the ideal gas laws, 
^d considering the uranium atoms to be doubly ionized. In mixing hydrogen 

hydrogen replace one of uranium and its two electrons 
to keep the pressure constant. 


5.8.1 


j^ff V8 Radius Ratio at Constants Vc^ume Fraction 


to fUa varying the ratio of the radius to which the fuel expands 

s >^adlus at a constant volume fraction, and hence fuel mass, was 

effect is shown graphically in Figure 5.35. This is the most 
important predominant parameter identified. 


5.8.2 


Keff VB Rftff at Constant R.R. and V.F. 


€u«i the relative effects of varying the regions of different 

r^aTnf constant fuel mass and radius ratio, a density-weighted effective 

raoius, Kg££ was defined 


R„ . -li-1 

S 2.rj^ At, 

This effect is shown graphically in Figure 5.36. 

5*8*3 Keff V8 Volume Fraction at Constant Radius Ratio 

4 ..4 “aturally intuitive feeling, that the more fuel the higher the 

criticality factor, was verified in one calculatlonal series. This result 
is shown in Figure 5.37. 
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5.8.4 Other Fnctors Influencing Kg |||: 

The ei:£ect of using U'”233 for fuel over U'”23S was shown to decrease 
critical mass by about a factor of two. Increasing the reflector thickness 
increases A coat'-benefit analysis will be necessary to obtain the 

optimum reflector thickness for an actual rocket engine. 

Diffusion vs Transport Theory Codes 

The initial calculations were performed with a single thermal group, 
no upscatter diffusion theory code, DISNEL(1®). Due to the uniform nature of 
the results of these calculations, showing no sensitivity to fuel arrangement 
in the cavity. The results were checked with a multi-thermal, full thermal 
upscatter S4 Transport Code SCAMP(ll) . There was a wide divergency of results 
indicating an Inadequacy in DISNBL for hot gas core work. Since transport 

computer time, a multi-thermal diffusion theory code 
MONA , using an energy format and cross section library identical to SCAMP 
WM used to see if the divergence was due to transport effects or upscatter 
effects in the thermal groups. MONA and SCAMP agreed to within l%Ak. 

Figure 5.38 shows a comparison of the three codes for a few models calculated. 

4 » trend of the calculations seem to lead to a confirmation 

or the feasibility of the gas core rocket engine concept. The results further 
point out the need for further study and optimisation of flowing gas cavity 
parameters for obtaining the best design for a nuclear rocket engine. 
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<400 


Fig. 5.38 Comparison of Computer Solutions — SCAMP and HQUA vs 




The fflnBibUity of tho nan corn roactor concept hnn boon strengthened 
by the experiments reported herein. It ennnot be stntod that absolute feoslblUty has 
been demonstrated, since the Interaction effects of thermodynamics (rodinnt and 
some convective heat transfer) with the fluid hydraulics hove not been export^ 
mentally studied except on a smnll scale (Reference TAPAwork), However, the 
nuclear coupling with the fluid flow patterns reported in the previous sections 
has been studied, both under room temporoturo conditions and at high tempera- 
ture. 


General conclusions that con be stated at this time are the following: 

1. The gas core continues to appear feasible. However, operation 
in the upflrlng (against acceleration or gravity) direction, as 
could be done for MHD applications but which would require thrust 
reversal for a rocket, makes the attclnlng of a critical configu- 
ration much easier. Protection of the cavity wells from burnout 
Is more difficult, however. 

2. If operation is to be in the downfiring direction, as with a 
rocket, then center gas injection in a direction other than 
downward from the top appears to offer significant enhancement of 
fuel density and hence of nuclear reactivity. 

3. A central "updraft” can be enhanced by shaping of the walls. Such 
^ updraft" along with alternate injection schemes as mentioned in 
(2) above also can result in significant enhancement of fuel density. 
These schemes are found essential when the center gas was much heavier 
than the outer gas. 

4. With dtandard center injection near the top, the theoretical^®^ 
predictions by Putre of density vs buoyancy number have been 
essentially verified experimentally. These do show that the 
large cavity, low thrust, high temperature operating design 
conditions will not be able to operate nuclearly unless 
schemes such as outlined in conclusions (2) and (3) above are 
employed. 

5. The hot hydrogen nuclear effect Is extremely significant and 
highly controlling. If anything, it would appear to add a 
high stability factor to the cavity operation. 

As to recommendations for areas of most critical Investigative need, the 
following seem most pertinent at this time: 

1. Start-up conditions with dust fuel appeared woefully lacking 
in the needed density of fuel when run at low pressure. High 
pressure testing needs to be done to determine if dust is a 
feasible start-up mechanism. 

2. The hot hydrogen effect was studied only somewhat cursorily, 
and found to be extremely significant. More attention and 
detail needs to be given to this unusual nuclear problem. 



The obflwved strong mixing at the edge pf the highly turbulent 
boundory layer leads to the conclusion that convection in this 
region may be extrmely significant. If indeed convection 

added to the radiative trnnafer model, these studies 
should be undertaken before further serious design work is done. 
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